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Abstract. We study a coupled deterministic system of vorticity evolution 
and salinity transport equations, with spatially correlated white noise on the 
boundary. This system may be considered as a model for gravity currents 
in oceanic fluids. The noise is due to uncertainty in salinity flux on fluid 
boundary. After transforming this system into a random dynamical system, 
we flrst obtain an asymptotic estimate of enstrophy evolution, and then show 
that the system is ergodic under suitable conditions on mean salinity input 
flux on the boundary, Prandtl number and covariance of the noise. 



1. Geophysical background 

A gravity current is the flow of one fluid within another driven by the gravitational 
force acting on the density difference between the fluids. Gravity currents occur in 
a wide variety of geophysical fluids. 

Oceanic gravity currents are of particular importance, as they are intimately related 
to the ocean's role in climate dynamics. The thermohalinc circulation in the ocean 
is strongly influenced by dense-water formation that takes place mainly in polar 
seas by cooling and in marginal seas by evaporation. Such dense water masses are 
released into the large-scale ocean circulation in the form of overflows, which are 
bottom gravity currents. 

We consider a two-dimensional model for oceanic gravity currents, in terms of the 
Navier-Stokes equation in vorticity form and the transport equation for salinity. 
The Neumann boundary conditions for this model involve a spatially correlated 
white noise due to uncertain salinity flux at the inlet boundary of the gravity 
currents. 

In the next section, we present the model and reformulate it as a random dynam- 
ical system, and then discuss the cocycle property and dissipativity of this model 
in §3 and §4, respectively. Main results on random attractors, enstrophy and er- 
godicity are in §5. Enstrophy is one half of the mean-square spatial integral of 
vorticity. Ergodicity implies that the time average for observables of the dynamical 
system approximates the statistical ensemble average, as long as the time interval 
is sufHciently long. 

2. Mathematical model 

Oceanic gravity currents are usually down a slope of small angle (order of a few 
degrees). We model the gravity currents in the downstream- vertical plane, and 
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ignore the variabihty in the cross-stream direction. This is an appropriate ap- 
proximate model for, e.g., the Red Sea overflow that flows along a long narrow 
channel that naturally restricts motion in the lateral planar plane 1 1 . In fact, we 
will ignore small slope angle and the rotation, both affect the following estimates 
non-essentially, i.e., our results below still hold with non-essential modification of 
constants in the estimates and in the conditions for the ergodicity. 
Thus we consider the gravity currents in the downstream- vertical (x, z)-plane. ft is 
composed of the Boussinesq equations for ocean fluid dynamics in terms of vorticity 
q{x,z,t), and the transport equation for oceanic salinity S{x,z,t) on the domain 
D = {{x,z) : < < 1}: 

qt + Jiq,i^) =Aq-Rad^S, 

(1) 1 

St + J{S,^)^—AS, 
Pr 

where 

q{x,z,t) = - All; 

is the vorticity in terms of stream function tJj, Pr is the Prandtl number and Ra 
is the Rayleigh number. Moreover, J(g, h) = gxh^ — gzh^ is the Jacobian operator 
and A = dxx + dzz is the Laplacian operator. All these equations are in non- 
dimensionalized forms. For the simplicity, we let Pr = 1. 

Note that the Laplacian operator A in the temperature and salinity transport 
equations is presumably dxx + with 5 being the aspect ratio, and kh, kv 

the horizontal and vertical diffusivities of salt, respectively. However, our energy- 
type estimates and the results below will not be essentially affected by taking a 
homogenized Laplacian operator A — dxx + dzz- AH our results would be true 
for this modified Laplacian. The effect of the rotation is parameterized in the 
magnitude of the viscosity and diffusivity terms as discussed in ^5] . 
The fluid boundary condition is no normal flow and free-slip on the whole boundary 

ip = Q, q = 0. 

The flux boundary conditions are assumed for the ocean salinity S. 
At the inlet boundary {a; = 0, 0<z<l} the flux is specified as: 

(2) dxS = F{z) + w{z,t), 

with F{z) being the mean freshwater flux, and the fluctuating part w{z, t) is usually 
of a shorter time scale than the response time scale of the oceanic mean salinity. 
So we neglect the autocorrelation time of this fluctuating process and thus assume 
that the noise is white in time. The spatially correlated white-in-time noise 'w{z,t) 
is described as the generalized time derivative of a Wiener process w{z,t) defined 
in a probability space (ri,F, P), with mean vector zero and covariance operator Q. 
On the outlet boundary {x = 1, < z < 1}: 

dxS = 0. 

At the top boundary 2; = f , and at the bottom boundary z = 0: 

dzS = Q. 

This is a system of deterministic partial differential equations with a stochastic 
boundary condition. 
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3. COCYCLE PROPERTY 

In this section we will show that 1^ has a unique solution, and by reformulating 
the model, we see it defines a cocycle or a random dynamical system. 
For the following we need some tools from the theory of partial differential equa- 
tions. 

Let W2{D) be the Sobolev space of functions on D with first generalized derivative 
in 1/2 (-D), the function space of square integrable functions on D with norm and 
inner product 



\W\\l2 = yj \u{x)\ dD j , {u,v)l2^ J u{x)v{x)dD, u, v £ L2{D). 
The space W2{D) is equipped with the norm 

\H\wi = \Ml2 + \\dxU\\L2 + \\dzU\\L2- 

o 

Motivated by the zero-boundary conditions of q we also introduce the space wl {D) 
which contains roughly speaking functions which are zero on the boundary dD of 
D. This space can be equipped with the norm 

(3) \\uUiiD) = \\dM\L2 + \\dA\L2- 

Similarly, we can define function spaces on the interval (0, 1) denoted by ^2(0, 1). 
Another Sobolev space is given by W2{D) which is a subspace of H^2^(D) consist- 
ing of functions h such that hdD = 0. A norm equivalent to the H/2^-norm on 
W2{D) is given by the right hand side of (PJ. For the subspace of functions in 
L2{D) having this property we will denote as L2{D). 

We reformulate the above stochastic initial-boundary value problem into a random 
dynamical system pp. For convenience, we introduce vector notation for unknown 
geophysical quantities 

(4) u^{q,S). 

Let w be a white noise in £2(0, 1) with finite trace of the covariance operator Q, 
and the Wiener process w{t) be defined on a probability space (ri,F, P). 

Now we choose an appropriate phase space H for this system. We assume that the 
mean salinity flux F G £2(0, 1). Note that 

^ / SdD ^ I + = constant. 

Jd Jo 

It is reasonable (see [S]) to assume that 

(5) / [F{z) + w{z,t)]dz^O, 







and thus SdD is constant in time and, without loss of generality, we may assume 
it is zero (otherwise, we subtract this constant from S): 



SdD = 0. 

D 



Thus S e L2{D), and we have the usual Poincare inequality for S. 
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Define the phase space 

H = L2{D) xL2{D). 
We rewrite the coupled system Q as: 
du 

(6) — + Au = Fi{u) + F2{u), u{0) = uoeH, 

at 

where 

Au 

Fi{u)[x, z] 

and 

F2{u)[x,z] 

The boundary and initial conditions are: 

0,on dD, 

0, ondD\{x = 0, < z < 1}, 
F{z) + w{z, t), on {a; = 0, < z < 1}, 

where n is the out unit normal vector of dD. The system © consists of determin- 
istic partial differential equations with stochastic Neumann boundary conditions. 



We now transform the above system (jSJ into a system of random partial differential 
equations (i.e., evolution equations with random coefficients) with homogeneous 
boundary conditions, whose solution map can be easily seen as a cocycle. Thus we 
can investigate this dynamics in the framework of random dynamical systems PP. 
Note that we have a non-homogenous stochastic boundary condition for salinity 5, 
so the first step is to homogenize this boundary condition. 

To this end, we need an Ornstein-Uhlenbeck stochastic process solving the linear 
differential equation 



~ ( -Ias ) ' 

-( -^('^'^^^ ]\x z] 



(7) 



q 
as 

dn 

as 

dx 



u{0) 



(8) 



aril . 
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with the fohowing same boundary conditions as for S, and zero initial condition 

9j,77i(i, 0,z,w) ^F{z) +w{z,t), 
d^mit, i,z,uj) =0, 
dzrii{t,x,0,uj) =0, 
(9) a,77i(t,a;,l,c^) =0, 

r]idD =0. 



Lemma 3.1. Suppose that the covariance Q has finite trace : tr^^Q < oo. Then 
the Ornstein-Uhlenbeck problem ^-^^ has a unique stationary solution in L2{D) 
generated by 

{t,uj) rii{9tuj). 

In fact, we can write down the sohition rji foUowing Da Prato and Zabczyck ji;-{lll2| 
as 



(10) 



T]i{t,x,z,uj) = (-A) f S{t- s)Af{X)ds 
Jo 



where / is the identity operator in L2{D), and TV is the solution operator to the 
elliptic boundary value problem Ah — A/i = with the boundary conditions for h 
the same as ryi, that is dh/dn = X on dD with hdD = 0, where n is the unit 
outer normal vector to dD and X is 

/ Fiz) + w{z,s) \ 






V 



Here A is chosen so that this elliptic boundary value problem has a unique solution. 
Since hdD = 0, we can choose A = 0. Moreover, S(t) is a strongly continuous 
semigroup, symbolically, e^*, that is, the generator of S{t) is A. 



(11) 



r]{t,x,z,uj) = 



Now we are ready to transform J^Jl into a random dynamical system in Hilbert 
space H. Define 



and recall 
Let 



(12) V :— u — r/. 

Then we obtain a random partial differential equation 



dn 

(13) — + Av^F,{v + r,{etu;j) + F2iv + rj{etLj)), v{0) = vo e H. 
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where 

Vx(t, 0, z, Lu), Vx{t, 1, z, w), Vz{t, X, 0, uj) and Vz{t, x, 1, lo) are now all zero vectors, i-e., 
homogenous boundary conditions, and the initial condition is still the same as for 
u 




However, because of the Jacobian, we will have to do nonlinear analysis on to 
resolve v. 



We introduce another space 



V = wl{D) X Wi{D) 



For sufhciently smooth functions v 
parts 



(q, S*), we can calculate via integration by 



(14) 




since now, we have only homogenous boundary conditions. 

Hence on the space V we can introduce a bilinear form a(-, •) which is continuous, 
symmetric and positive 




This bilinear form defines a unique linear continuous operator A : V 
that {Au,v) — a{u,v). 

Recall 



V such 




and 




Lemma 3.2. The operator Fi : V 



H is continuous. In particular, we have 



(Fi(w),u> =0. 



Proof. We have a constant ci > such that 



(15) 



IIV'llwl(D) < ClMwiiD) 
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for any q £ W2{D) which follows straight forwardly by regularity properties of a 
linear elliptic boundary problem. Note that W2 is a Sobolev space with respect to 
the third derivatives. Hence we get: 

||^(5',V')I|l2 < sup {\d.^tl){x,z)\^\d^'4){x,z)\)-x 

X (^J \d^S{x,z)\ + \d,S{x,z)\dD 
The second factor on the right hand side is bounded by 



d^S{x,z)\^dDj +\^J \d,Six,z)\'dDj <\\u\\v. 

On account of the Sobolev embedding Lemma, we have some positive constants 
C2, C3 such that 

sup {\d^ip{x,z)\ + \d^ip{x,z)\) < czllVV'llw/KD) < cahWwiiD) <C3\\u\\v. 

(x.z)£D 

Hence we have a positive constant C4 such that 

\\J{S,n\L,<c4u\\l 

for u . 

We now show that 

(J(5,V),5) =0. 
We obtain via integration by parts 

d^Sd.ipSdD^ [ d.Sd^ipSdD 

JD 

= - I d^,S^/jSdD+ [ dl^S^pSdD- [ d^Si^d^S dD + [ d^Sipd^S dD 

JD JD JD JD 

+ [ d^S^S\lzldx- [ d^Si:S\%zldz = 

"'(0,1) "'(0,1) 

because ip is zero on the boundary dD. This relation is true for a set of sufficiently 
smooth functions ip, S which are dense in wli^) x W2{D). By^the continuity of i^i, 
as just shown in Lemma IT^ we can extend this property to W^iD) x W2{D). □ 



Lemma 3.3. The following estimate holds 

\\F2{u)\\l2 < C5\\u\\v. 

for some positive constant C5. 

Proof. By simple calculation, the proof is obtained. □ 

We have obtained a differential equation without white noise but with random 
coefficients. Such a differential equation can be treated sample-wise for any sample 
Lu. We are looking for solutions in 

veC{[0,T];H)nL^{0,T;V), 

for all r > 0. If we can solve this equation then u := v + r] defines a solution version 
of ®. 



yENA PEARL BONGOLAN- WALSH, JINQIAO DUAN, HONGJUN GAO, TAMAY OZGOKMEN, PAUL FISCHER, AND TRAIAN ILIESC 

For the weU posedness of the problem we now have the following result. 

Theorem 3.4. (Well- Posedness j For any time t > 0, there exists a unique 
solution of in C{[0,t]; H) 0X2(0, t;^). In particular, the solution mapping 

R+ xQx H 3 {t,uj,vo) v{t) £ H 

is measurable in its arguments and the solution mapping H 3 vq ^ v(t) E H is 
continuous. 

Proof. By the properties of A and Fi (see Lemma \3.2\i . the random differential 
equation H13() is essentially similar to the 2 dimensional Navier Stokes equation. 
Note that F2 is only an afhne mapping. Hence we have existence and uniqueness 
and the above regularity assertions. □ 

On account of the transformation (|12|l . we find that (O also has a unique solution. 
Since the solution mapping 

IR+ X Q X H 3 {t^LO, Vq) — > v{t, LU, Vq) —: ip{t, lu, Vq) G H 

is well defined, we can introduce a random dynamical system. On we can define 
a shift operator 9t on the paths of the Wiener process that pushes our noise: 

w{-,9tuj) = w{- + t,uj) - w{t,uj) for t e R 

which is called the Wiener shift. Then {^tjteffi forms a flow which is ergodic for the 
probability measure P. The properties of the solution mapping cause the following 
relations 

(p(t + T, U!, u) = (p(t, OrUJ, ip{T, U!, u)) for t, T > 
</?(0, U!,u) — U 

for any u £ and u € H . This property is called the cocycle property of (p which 
is important to study the dynamics of random systems. It is a generalization of 
the semigroup property. The cocycle if together with the flow 6 forms a random 
dynamical system. 

4. DiSSIPATIVITY 

In this section we show that the random dynamical system H13|) for gravity currents 
is dissipative, in the sense that it has an absorbing (random) set. This means that 
the solution v is contained in a particular region of the phase space H after a suf- 
ficiently long time. This dissipativity will help us to obtain asymptotic estimates 
of the enstrophy and salinity evolution. Dynamical properties that follow from this 
dissipativity will be considered in the next section. In particular, we will show that 
the system has a random attractor, and is ergodic under suitable conditions. 

We introduce the spaces 

H = L2{D) 
V^WliD). 

We also choose a subset of dynamical variables of our system (Q. 
(16) v^S-r^i. 
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To calculate the energy inequality for v, we apply the chain rule to HwHl/- We 
obtain by Lemma 18.21 

(17) llNl| + 2||V.||L 

-2(J(77i,V),S). 

The expression Vw is defined by (S7x,zv)- We now can estimate the term on the 
right hand side. 

By the Cauchy inequality, integration by parts and Poincare inequality Ai||(7||l2 < 
||Vg||L2 for q e Wl{D) and X2ML2 < W^vh^ for v G WUD), we have 

(18) |ll^ll| + l|V^||L<A?||Vg|nhif. 
For q, we have the following estimate 

(19) |lklP + l|Vgf <i?a2(A2||V^f + hif). 
From (^HJ and we have 



Definition 4.1. A random set B = {B(uj)}uj^n consisting of closed hounded sets 
B{uj) is called absorbing for a random dynamical system Lp if we have for any 
random set D — {D(a;)}(jgo, D{uj) G H bounded, such that t — > supj^g^j^g^j^-j \\y\\H 
has a subexponential growth for t — s- ±00 

(p{t,uj,D{Lu)) C B{9tLj) fort>to{D,Lu) 

^it,e-tUJ,D{e^tio)) CB{uj) fort>toiD,Lu). 

B is called forward invariant if 

(p{t,u),uo) £ B{6tiLj) if uq £ B{lo) fort>0. 

Although u is not a random dynamical system in the strong sense we can also show 
dissipativity in the sense of the above definition. 

Lemma 4.2. Let ip{t,uj,vo) G H for vq £ H be defined in and 

^ 2A?]E||r;i|2 >0. 



Ra^Xl 

Then the closed ball B{0, Ri{uj)) with radius 

T 

A 



" 1 



i?iM = 2/ e^^^WmWdr 



is forward invariant and absorbing. 

The proof of this lemma follows by integration of (|2U|I . 

For the applications in the next section we need that the elements which are con- 
tained in the absorbing set satisfy a particular regularity. To this end we introduce 
the function space 

:= {u e H : Ml := \\A^u\\j,<oo} 
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where s G M. The operator A'^ is the s-th power of the positive and symmetric op- 
erator A. Note that these spaces are embedded in the Slobodeckij spaces H'^, s > 0. 
The norm of these spaces is denoted by || • \\h^- This norm can be found in Egorov 
and Shubin 7 , Page 118. But we do not need this norm expUcitly. We only men- 
tion that on the norm || ■ ||s of is equivalent to the norm of W for < s, see 

Our goal is it to show that v{l, us, D) is a bounded set in TL'^ for some s > 0. This 
property causes the complete continuity of the mapping v{l,uj, ■). We now derive 
a differential inequality for t||ti(i)||^. By the chain rule we have 

Note that for the embedding constant Cg between Ti^ and V 



f\\v\\lds<cl f\ 
JO Jo 



WvWyds for s < 1 

such that the left hand side is bounded if the initial conditions vq are contained in 
a bounded set in H. The second term in the above formula can be expressed as 
followed: 

t^^{Aiv,Aiv)H =2t{j^v,A'v)H - -2t{Au,A'v)H + 2t{Fiiv + 7^ietiu)),A'v)H 
+ 2t{F2{v + Tj{0tiu)),A'v)H- 

We have 

iAv,A%)H ^ \\a'^+^v\\h = \\v\\i+,. 

Similar to the argument of |21| and the estimate for the existence of absorbing, and 
applying some embedding theorems, see Temam Page 12 we have got 

\\v\\l<Cit,\\vo\\H, sup fort e [0,1]. 

te[o,i] 

By the results of [H] and PT]. we know 

sup ||r?i||i>(A=) < C{trL2Q) < oo, for < s < ^. 
te[o,i] 4 

By now, our estimates allow us to write down the main assertion with respect to 
the dissipativity of this section. 

Theorem 4.3. For the random dynamical system generated by ^Ici}} . there exists a 
compact random set B = {B{uj)}^^f,-i which satisfies Definition \4.1\ 

We define 

(22) B{uj) = (p{l,e^iuj,B{0, R{9 ^luj))) C H', < s < i. 

In particular, Ti" is compactly embedded in H. 



5. Random dynamics: Enstrophy and ergodicity 

In this section we will apply the dissipativity result of the last section to analyse 
the dynamical behavior of the random dynamical system However, it will be 
enough to analyse the transformed random dynamical system generated by (|13|l . 
By the transformation (fT^ . we have all these qualitative properties to the original 
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gravity current system l(T]l. 

We will consider the following dynamical behavior: random attractors, asymptotic 
evolution of enstrophy and mean-square norm of salinity profile, and ergodicity. 
Enstrophy is defined as one half of the mean-square integral of vorticity. Ergodicity 
implies that the time average for observables of the dynamical system approximates 
the statistical ensemble average, as long as the time interval is sufficiently long. 
We first consider random attractors. We recall the following basic concept; see, for 
instance, Flandoli and Schmalfufi 

Definition 5.1. Letip be a random dynamical system. A random set A = {A{uj)}^(^q 
consisting of compact nonempty sets A{uj) is called random global attractor if for 
any random bounded set D we have for the limit in probability 

(P) lim distH{^p{t,u;,D{Lj)),A{etLj)) = 

t — ^oo 

and 

(p{t,uj,A{uj)) ^ A{etuj) 

any t > and w G 51. 

The essential long-time behavior of a random system is captured by a random 
attractor. In the last section we showed that the dynamical system (p generated by 
is dissipative which means that there exists a random set B satisfying H21() . In 
addition, this set is compact. We now recall and adapt the following theorem from 

m 

Theorem 5.2. Let ip be a random dynamical dynamical system on the state space H 
which is a separable Banach space such that x ip{t,oj,x) is continuous. Suppose 
that B is a set ensuring the dissipativity given in Definition In addition, B 
has a subexponential growth (see Definition \4.1\j and is regular (compact). Then 
the dynamical system ip has a random attractor. 

This theorem can be applied to the random dynamical system ip generated by the 
stochastic differential equation itT^ . Indeed, all the assumptions are satisfied. The 
set B is defined in Theorem 14.31 Its subexponential growth follows from B{uj) C 
B{0, R{u!)) where the radius R{u!) has been introduced in the last section. Note 
that is a continuous random dynamical system; see Theorem ll'i.41 Thus ip has 
a random attractor. By the transformation (|12|l . this is also true for the original 
gravity currents system 0. 

Corollary 5.3. fRandom Attractor^ The gravity currents system {^j has a 
random attractor. 

Now we consider random fixed point (stationary state) and ergodicity. 

Definition 5.4. A random variable v* : CI ^ H is defined to be a random fixed 
point for a random dynamical system if 

(p{t,uj,v*{uj)) = V*{0tUj) 

for t > and lo £ fl. A random fixed point v* is called exponentially attracting if 

lim \\ip{t,LU,x)-v*{etLo)\\H = 

for any x ^ H and cu ^ il. 
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Sufficient conditions for the existence of random fixed points are given in Schmalfufi 
|14j . We here formulate a simpler version of this theorem and it is appropriate for 
our system here. 

Theorem 5.5. (Handom Fixed Point Theorem^ Let f he a random dynamical 
system and suppose that B is a forward invariant complete set. In addition, B has 
a subexponential growth, see Definition \4.1\ Suppose that the following contraction 
conditions holds: 



where the expectation o/logfc denoted by Elogfc < 0. Then ip has a unique random 
fixed point in B which is exponentially attracting. 

This theorem can be considered as a random version of the Banach fixed point 
theorem. The contraction condition is formulated in the mean for the right hand 
side of 

Theorem 5.6. (TJnique Random Stationary Statej Assume that the salinity 
boundary flux data \\F\\l2, the Prandtl number Pr, and the trace of the covariance 
for the noise tr^^Q are sufficiently small. Then the random dynamical system 
generated by |^ has a unique exponentially attracting random stationary state. 

Note that if we take into account the horizontal and vertical diffusivities kh^i^v 
in the Laplacian operator, the above "smallness" condition needs to be slightly 
modified. 

Here we only give a short sketch of the proof. Let us suppose for a while that B is 
given by the ball i?(0, R) introduced in Lemma ??. Suppose that the data in the 
assumption of the lemma are small and v is large. Then it follows that Ei? is also 
small. To calculate the contraction condition we have to calculate w, Vi{uj)) — 
93(1, uj, V2{'.^))\W for arbitrary random variables vi, V2 £ B. By Lemma lH.2l we have 
that 

(J(gi,Vi) - J(gi,'0i),gi -92) < C22|ki - q2\\wi + C23lki||^i(£,)||'7i - q2\\l^ 

where the constant C23 can be chosen sufficiently small if v is large. On account of 
the fact that the other expressions allow similar estimates and that F2 is linear we 
obtain: 

^\\f{t,UJ,Vi{uj)) - (p{t,UJ,V2{uj))\\% 

< {-a' + C23\\v{t,uj,V2{uj))\\v)\\(p{t,u;,vi{uj)) ~ if{t,uj,V2{uj))\\H 

for some positive a' depending on ly. From this inequality and the Gronwall lemma 
it follows that the contraction condition H23|) is satisfied if 



(23) 



sup 



\ip{l,uj,vi) - kp{\,ij:,V2)\H 



< k{uj) 




But by the energy inequality this property is satisfied if the Ei? and E||z||y is suf- 
ficiently small which follows from the assumptions. 
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Let now B be the random set defined in H22|l . Since the set B introduced in H22|l 
is absorbing any state the fixed point v* is contained in this B. In addition v* 
attracts any state from H and not only states from B. 

The uniqueness of this random fixed point impHes ergodicity. We will comment on 
this issue at the end of this section. 



By the well-posedness Theorem ll-i.4l we know that the stochastic evoltuion equation 
© has a unique solution. The solution is a Markov process. We can define the 
associated Markov operators T{t) for i > 0, as discussed in E]. Moreover, 
{T{t)}t>o forms a semigroup. 

Let be the set of probability distributions /i with finite energy, i.e.. 



H 



Then the distribution of the solution u{t) (at time t) of the stochastic evolution 
equation © is given by 

where the distribution of the initial data is contained in M^. 

We note that the expectation of the solution can be expressed in terms of 

this distribution T(t)^iQ: 



E\\u{t)\\i,= / \\u\\Wr{t)^io. 

Jh 

We can derive the following energy inequality in the mean, using our earlier esti- 
mates. 

Corollary 5.7. The dynamical quantity u = (q, S) of the coupled atmospheric- 
ocean system satisfy the estimate 

E\\u{t)\\l + aE f \\uiT)\\ldT <E\\uo\\jj + tc24 + ttTL,Q, 
Jo 

where the positive constants C24 and a depend on physical data F{x), Pr and Ra. 

By the Gronwall inequality, we further obtain the following result about the asymp- 
totic mean-square estimate. Recall that Jjj q^dD is the enstrophy for gravity 
currents. The following theorem gives asymptotic estimate for both enstrophy and 
mean-square salinity E S^dD. 

Theorem 5.8. (Enstrophy and mean-square salinityj For the expectation of 
the dynamical quantity u = {q, S) of the coupled atmospheric- ocean system, we have 
the asymptotic estimate 

limsupE||M(i)|||^ = limsup f \\u\\'jfdT{t)^Q < "^^"^ ^''^2*3 

t— >CX3 t—too J H ^'^5 

if the initial distribution fiQ of the random initial condition Uoit-^) is contained in 
M"^ . Here C25 > depends only on physical data. 

By the estimates of Corollary 15.71 we are able to use the well known Krylov- 
Bogolyubov procedure to conclude the existence of invariant measures of the Markov 
semigroup T{t). 
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Corollary 5.9. The semigroup of Markov operators {T{t)}t>o possesses an invari- 
ant distribution fjLi in M"^ : 

T{t)fii = fii for t >0. 

In fact, the limit points of 

(i /V(T)M0dTj 

L ^ Jo J t>o 

for t — » 00 are invariant distributions. The existence of such limit points follows 
from the estimate in Corollary 15. 71 

In some situations, the invariant measure may be unique. For example, the unique 
random fixed point in Theorem 15.61 is defined by a random variable u*{uj) = 
v*{ll!) + rjiuj). This random variable corresponds to a unique invariant measure 
of the Markov semigroup T{t). More specifically, this unique invariant measure is 
the expectation of the Dirac measure with the random variable as the random mass 
point 

Because the uniqueness of invariant measure implies ergodicity jli-il . we conclude 
that the gravity currents model Q is ergodic under the suitable conditions in 
Theorem l5.6l for physical data and random noise. We reformulate it as the following 
ergodicity principle. 

Theorem 5.10. (^rgodicityj Assume that the salinity boundary flux data \\F\\l2, 
the Prandtl number Pr, and the trace of the covariance for the noise tr^^Q are 
sufficiently small. Then the gravity currents system 0) is ergodic, namely, for 
any observable of the gravity currents, its time average approximates the statistical 
ensemble average, as long as the time interval is sufficiently long. 

In the regime of ergodicity, the gravity currents system can be numerically simulated 
with (almost surely) one random sample. This is what we call ergodicity-based 
numerical simulation of stochastic systems, in contrast to sample- wise ("Monte 
Carlo") simulations. 
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